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We study the dissipation in a system consisting of a small metallic island coupled to a gate 
electrode and to a massive reservoir via single tunneling junction. The dissipation of energy is 
caused by a slowly osciUating gate voltage. We compute it in the regimes of weak and strong 
Coulomb blockade. We focus on the regime of not very low temperatures when electron coherence 
can be neglected but quantum fluctuations of charge are strong due to Coulomb interaction. The 
answers assume a particularly transparent form while expressed in terms of specially chosen physical 
observables. We discovered that the dissipation rate is given by a universal expression in both 
limiting cases. 
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I. INTRODUCTION 



The phenomenon of Coulomb blockade has become an 
excellent tool for observation of interaction effects in sin- 
gle electron devices. Theoretical means for its explo- 
ration are well developed and versatile J^i^i^i^i^ The sim- 
plest mesoscopic system displaying Coulomb blockade is 
a single electron box (SEB). The properties of such a 
system are essentially affected by electron coherence and 
interaction. Our work is motivated by a considerable re- 
cent theoretical and experimental interest in the relation 
between dissipation and resistance of this device in vari- 
ous parametric regimes>^i^i^ i^°i^^i^^ 

The set-up is as follows (see Figll}. Metallic island is 
coupled to an equilibrium electron reservoir via tunneling 
junction. The island is also coupled capacitively to the 
gate electrode. The potential of the island is controlled 
by the voltage Ug of the gate electrode. The physics 
of the system is governed by several energy scales: the 
Thouless energy of an island -Exh, the charging energy 
Ec, and the mean level spacing S. Throughout the paper 
the Thouless energy is considered to be the largest scale 
in the problem. This allows us to treat the metallic is- 
land as a zero dimensional object with vanishing internal 
resistance. The dimcnsionless conductance of a tunneling 
junction g is an additional control parameter. 

Initially, the main quantity of interest in a Coulomb 
blockaded SEB was its effective capacitance: dQ/dUg, 
where Q is the average charge of a islandj ^'^i^^'^^i^^'^^i^^ 
Paper [7| however sparked both theoretical and exper- 
imental attention to the dynamic response functions of 
such a set-up<^i^ii2iiiiii^'^i2i It is worthwhile to mention 
that the system docs not allow for conductance measure- 
ments since there is no DC-transport. This way an es- 
sential dynamic characteristic becomes the set-up admit- 
tance, which is a current response to an AC-gate voltage 
Ug{t) = Uq + Uui cosujt. As it is well-known, the real part 
of admittance determines energy dissipation in an elec- 
tric circuit. Classically, the average energy dissipation 




FIG. 1: Measurement of resistance Rq. The SEB is subjected 
to a constant gate voltage Uo The dissipative current through 
the tunneling contact is caused by a weak AC voltage U{t). 



rate of a single electron box is given as follows 



u^C^gR]U^\^, 



R 



hui-^gE^, (1) 



where Cg denotes the gate capacitance, e - the electron 
charge, and h = 2-111% - the Planck constant. Expres- 
sion ([T]) presents us with a natural way of extracting the 
resistance of a system from its dissipation power. The re- 
sistance of a classical system is thus fully determined by 
the tunneling conductance of the contact via Kirchhoff's 
law: R = h/ {e^g). The question one asks is how quantum 
effects such as electron coherence and interaction change 
this result? One expects that correct quantum dissipa- 
tion is going to give generalized quantum resistance. The 
obvious stumbling block one foresees is that only combi- 
nation of two observables: C^R can be extracted from the 
dissipation power rather than just R. For the case of fully 
coherent SEB this key difficulty was resolved in Ref. [31 ■ 
It was shown that the energy dissipation rate can be 
factorized in accordance with its classical appearance ([1]) 
but the definition of physical quantities comprising it be- 
comes different. Geometrical capacitance Cg should be 
substituted by a new observable: mesoscopic capacitance 
Cp. This leads to the establishment of another observ- 
able: charge relaxation resistance i?„ such that R —f Rq 



in Eq. ((T|). Charge relaxation resistance of a coherent sys- 
tem differs drastically from its classical counterpart. In 
particular, as shown in Ref. 0], the charge relaxation re- 
sistance of a single channel junction doesn't depend on its 
transmission. The admittance in the quasi-static regime 
was investigated in the recent experiment by Gabelli et 
aliii The measurements were performed at low tempera- 
tures T < S when the system could be regarded as coher- 
ent. The question that has remained unattended by the 
theory is what happens to dissipation and resistance at 
transient temperatures when thermal fluctuations smear 
out electron coherence but electron-electron interaction 
is strong? The recent experiment by Persson et ali^ ex- 
plored the energy dissipation rate at these transient tem- 
peratures. 



Motivated by the experiment [13 we address the same 
question from the theoretical point of view. We study 
the energy dissipation rate of a single electron box in 
the so-called 'interactions without coherence' regime. It 
corresponds to the following hierarchy of energy scales: 
i?Th 3> i?c > r max{i5, g(5}. This temperature 
regime is such that keeps electrons strongly correlated 
(T <C Ec), yet allows to discard electron coherence 
(T 3> max{(5, (7(5})i22iS We compute the energy dissipa- 
tion rate and the SEB admittance in the limits of large 

3> 1) and small (.g <C 1) dimensionless tunneling con- 
ductance of the jimction. 

We consider a multichannel junction but the conduc- 
tance of each channel is assumed to be small gch ^ 1- 
Then, the physics of the system is most adequately de- 
scribed in the framework of Ambegaokar-Eckern-Schon 
(AES) effective actioui^ Our results lead to the gener- 
alization of classical result H]). We found that at w ^ 
the average energy dissipation rate can be factorized in 
both g ^ I and <? ^ 1 limits as 



h 



eV(r)' 



(2) 



in complete analogy with classical expression ([T]). Here, 
Rq{T) and Cg{T) are identified as charge relaxation re- 
sistance and rcnormalized gate capacitance, respectively. 
It is worthwhile to mention that the physical observables 
g'(T) and Cg{T) are defined universally for any value of 
dimensionless conductance g. It allows us to suggest that 
Eq. ([2]) remains valid for arbitrary value of g. 

In order to explain physics behind quantities g'(T) and 
Cg{T), it is useful to consider a single electron transis- 
tor (SET) rather than SEB (see Fig.[2|). In the absence 
of DC-voltage between left and right reservoirs a SET 
represents essentially a SEB except different definition of 
the parameter g. Then. g'{T) is the very quantity that 
determines the SET conductance. The rcnormalized gate 
capacitance Cg (T) is very different from the effective ca- 
pacitance dQ/dUo. In fact, Cg{T) = dq'{T)/dUo, where 
q'{T) is the physical observable introduced recently in 
Ref. [2^ to describe the 0-angle renormalization in the 
Coulomb blockade problem. The quantity q' is deter- 
mined not only by the average charge Q but also by the 
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FIG. 2: Measurement of conductance. The SET is subjected 
to a constant gate voltage Ug and constant bias U . 



anti-symmetrized (so-called, quantum) current noise in a 
SET. 

The paper is organized as follows. Section II is used to 
introduce AES model. Sections III and IV are devoted 
to dissipation in the weak [g ^ 1) and strong [g 1) 
coupling regimes. Section V is devoted to discussion. 



II. FORMALISM 

A. Hamiltonian 

A single electron box is described by the Hamiltonian 

H = Ho+H, + Ht, (3) 

where Hq describes free electrons in the lead and the 
island. He describes Coulomb interaction of carriers in 
the island, and Ht describes the tunneling. 



(4) 



Here, operators (dj,) create a carrier in the lead (is- 
land). 



Ht^^ tkaa\da + h.c. 



(5) 



The charging Hamiltonian of electrons in the box is taken 
in the capacitivc form: 



d-q) 



(6) 



Here, Ec = e^/(2C) denotes the charging energy, and 
q — CgUg/e the gate charge, fid is an operator of a 
particle number in the island: 



nd 



(7) 



It is convenient to introduce hermitean matrices: 



gkk' = (27r)2 



|„„, = (2^)2 5{e(^me'-'))\'''Y.^iA^k^)tko.', (9) 
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the first of them acting in the Hilbert space of the states 
of the lead, the second - in the space of the islands 
states. The energies e^''^ e^'') are accounted for with re- 
spect to the Fermi level, and the delta-functions should 
be smoothed on the scale 6E, such that 6 <C SE ^ T. 

The eigenstates of g (g) describe the 'channel states' 
in the lead (island), while the transmittances of the cor- 
responding channels Ty are related to the eigenvalues g-y. 
Note, that in general the rank of the matrix g differs from 
that of the matrix g, so that the numbers of eigenvalues 
are also different. This difference is, however, irrelevant, 
since it stems from the "closed channels" with w 0, 
i.e., the states strongly locahzed either within the lead, 
or within the island. The effective 'channel conductance' 
^ch and the effective number of open channels A^ch can 
be defined aa^ 



5ch 



tr(g') 
trg 



ch 



(trg) 
tr(g2) 



(10) 



In general case the effective action can be written as a 
sum of terms, proportional to tr(g'^), over all integer k 
(see ref.jl^). The problem is considerably simplified in 
the tunnel case, when 



5ch < 1- 



(11) 



In the present paper we assume this condition to be sat- 
isfied. Then all terms with fc > 1 can be neglected, com- 
pared to the leading term with k = 1 and the standard 
form of the AES action can be easily reproduced. In 
particular, the classical dimensionless conductance of the 
junction is expressed as: 



g = trg = trg = .gch^ch- 



(12) 



Note, that under the condition pT|) g still can be large, 
if the number of channels A'ch ^ 1 is sufficiently large. 

Throughout the paper we keep the units such that h = 
e = 1, except for the final results. 



B. Conductance and dissipation 

To study the electric properties of a system we compute 
energy dissipation caused by slow oscillations of external 
gate voltage Ug{t) = Uo + Ui^ coaojt. 

The average energy dissipation rate can be found fol- 
lowing the standard schemei^l 



dE / '^^ \ ^'-'g 
~ ~dt ^ \6U^/~dr 



SH \ dUn 



(13) 



Here, E is the energy of the system, H is given by ([3]) and 
angular brackets denote full quantum statistical average. 
Since 



/ 



SH 



Co 

—S-U 

c ' 



(14) 



the energy dissipation is determined by a response of the 
electron density in the island to the time-dependent gate 
voltage Ug{t). Therefore, it can be found via Callen- 
Welton fluctuation-dissipation theorem)^ 



(J2 

= ^wlmn«(c^)|i7^|2. 



(15) 



Here, n^(cL') is the retarded electron polarization opera- 
tor: 

H^(t) =^e(^)([nd(^),nd(0)]), rid = E4d„ (16) 



with 6(t) denoting Hcaviside step function. 

We are interested in the quasi-static regime u; — > 0. 
Then, as it will be proven below, the polarization op- 
erator n^(a;) is possible to expand in regular series in 



n^(w) = ^o(T) + iujni{T) + e)(tj2), 



(17) 



where both ttq (T) and tti (T) are real functions of tem- 
perature and other SEB parameters. Then the energy 
dissipation rate is solely determined by the linear coeffi- 
cient 7ri(T) and acquires Ohmic form: 



W^ = —A{T)\U^\', A{T) = -^MT). (18) 

The SEB admittance g{Lo) which is the linear reponse 
of an AC-current to AC-gate voltage U^^: G{(^) = 
Ii^/Ui^, is related to the polarization operator (see Ap- 
pendix): 



g{oj) = -iLuCg{i + n^{uj)/c). 



(19) 



As expected, the energy dissipation rate is proportional 
to the real part of the admittance: Wui ~ Keg{uj). The 
static part of the polarization operator 11^ (uj) is deter- 
mined by the effective capacitance dQ/dUo as 



C dQ 

'c'gdUo 



(20) 



where Q = {fid) denotes the average charge on the island. 
We mention that Eq. (|20p is analogous to the well-known 
Ward identity which relates static polarization operator 
and compressibility)^ Using Eqs. pT|) - (P0|) . we can es- 
tablish the following result: 

Giu;) = -zc.^ + §-AiT)u;^ + 0{lo^) (21) 

which is a quantum generalization of the classical relation 

g{uj) = ~iujCg + CgCRijj'^ + 0{uj^). (22) 

Therefore, both the admittance and the energy dissipa- 
tion rate are determined by the polarization operator 
H''^(a;) which involves one imknown function tti{T) in 
the quasi-static regime. 
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C. AES model 

The condition (fTTj) vahdates the use of AES-effective 
actional which describes the physics of the set-up in 
terms of a single quantum phase <p(t) fluctuating in Mat- 
subara time r: 

Saes = Sd + Sg + Sc- (23) 

Here, Sd is the dissipative part of the action in the stan- 
dard form: 

Sd = -^ /''a(ri2)e*'^(^^)-^'^(^^) dr.dr^, 

(24) 



a{T) 



Jo 

J.2 



T 



sin^ ttTt 



-E 



OJn e 



where (3 = l/T, T12 = ti — T2, ujn = ^irTn and g is 
defined by (|12p and stands for the dimensionless (in units 
e^/h) conductance of the tunnel junction. The term Sg 
represents a coupling with the gate voltage Uq: 



Sn 



-iq I (pdr 



-2TrqWi. 



(25) 



Here, integer W is the winding number of a field (pir) 
which appears through the constraint 



ifiP) - ifiO) = 2ttW. 



(26) 



Non-zero value of Sg appears for topologically non-trivial 
field configurations only. The charging part of the action 
is as follows 

r0 



Sc 



(27) 



Physically, time derivative of the phase variable ip de- 
scribes voltage fiuctuations in a SEB. We emphasize that 
AES-action is valid for any value of g. We work in 
the regime T <^ Ec- Charging term Sc is thus always 
small providing a natural ultraviolet cut-off in the the- 
ory: A = gEc. 

Our aim is to compute the polarization operator (jl6[) 
which, according to Eqs ([T8|) and (fT9|) . determines the 
energy dissipation and admittance. Therefore, we need 
to express initial observables cast in terms of fermionic 
operators through correlators of bosonic field (^(r). This 
is done in part A of the appendix by employing Keldysh 
formalism. The polarization operator Il^{ui) then can be 
obtained by analytical continuation iujn —^co + iOoi the 
following phase correlator in Matsubara basis: 



(28) 



Here %- denotes time ordering. So far we made no as- 
sumptions about the value of p. The AES model is how- 
ever impossible to tackle for arbitrary gs due to highly 
non-linear form of the dissipative term. In the next chap- 
ter of the paper we restrict our attention to the case of 
large dimensionless conductance 5 3> 1; the quantity 1/g 
then becomes an expansion parameter of perturbation 
theory. 



III. WEAK COUPLING REGIME, g > 1 

A. Perturbation theory 

To expand the polarization operator Il{iuj„) in powers 
of l/g it is convenient to use the Matsubara frequency 
representation. 



(29) 



Then, the quadratic part of AES action assumes the 
form: 



4'L=5E(-+^«1|^"^ (30) 



It determines the propagator of the ip field as 



1 



g |7i| + 27:^Tny{gEc) 



(31) 



Evaluation of the polarization operator at the tree level 
yields 



njiUJn) _ _2n\lJn\ , 2^ 

C2 " g 



(32) 



Performing standard one-loop calculations one finds 

' 2, 
1 + - In 



H(..„)_ 2.K|^^ , 2,„.gi^.e^-,^^(^2). (33) 



C2 



g 2'K^T 



With the help of the rcnormalization group analysis this 
result can be written aa^ 

'^ = '-'-^^Oiul). (34) 
Here, g{T) is given by 



g(T) g - 2 n^^-A^ 



(35) 



with 7 « 0.577 being Euler's constant. Eq. ([55]) describes 
the well-known one-loop temperature rcnormalization of 
the coupling constant 1^ 



B. Instantons 

So far the phenomenon of Coulomb blockade i.e. de- 
pendence on q, is completely absent in all our expres- 
sions for polarization operator. To catch it we have to 
take into account instanton solutions of AES actioui^i^^ 
Korshunov's instantons read 



Jipw{r\{z^}) _ 



\W\ 

n 



1 - z*e 



*„2mTT 



sgn W 



(36) 
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Here, Za is a set of arbitrary complex numbers. Positive 
values of winding numbers W are assigned to instantons 
with 1 0a I < 1 and negative ones to anti- instantons with 
1 0a I > 1- On the classical solutions ([55]) the dissipative 
Sd and topological Sg part of AES-action becomes 

.9: 



SdVPw] + Sg[ipw] 



\W\ - 2-KWqi. 



(37) 



It is finite and independent of ZaS. These parameters are 
zero-modes. The charging term though does depend on 
them; 

Sc['fw]^ -j^l^ , • (38) 

a,o " 

Thus ZaS can only be viewed as approximate zero modes 
and the instanton configurations with \za\ — *■ 1 are su- 
pressed. 

As it is clear from Eq. ([57]) every instanton brings a 
small factor e~^/^ to any observable we want to compute. 
In what follows, we restrict ourselves to one-instanton 
{W = ±1) contribution only. 



C. Instanton correction to the polarization 
operator 

To get the instanton contribution to the polarization 
operator we need to compute one-instanton correction 
to the correlator (7^(p(t)(^(0)). Up to the one-instanton 
contributions we find 



n(ia;n) 

Q2 



W=±l 

I + 11, 



W=±l 



(39) 



(40) 



where {tp^p)^^ = Jq''((p(t)(^(0)) exp(iu)„T)(iT and 
Jw 

^^^)iw) ^ 1 /■ p{T)meM-SAEs]. 
^0 Jw 

Here, the subscript W at the integral sign means that 
functional integration is performed over phase configu- 
rations obeying the boundary condition (|26p . The first 
term I in Eq. ((551) represents the renormalization of the 
partition function due to instantons. The second term 
II is the contribution of the instanton solutions ip±i into 
the correlation function itself. The rcnormalizcd parti- 
tion function read a-'^^'^^i^^'^^ 



1 



E 

W=±l 



~Z^ 



cos27rg. (41) 



The contribution II consists of two terms: 

H= > ■ (w)L!r'= E ^^w^w)Z^ 

W=±l 



E 

W=±l 



(42) 



where the first term is a correlator of classical field con- 
figurations (|36p averaged over zero modes Za and the sec- 
ond term comes from fluctuations of phase p around the 
classical solution ipw- As shown in Appendix B the lat- 
ter term in (|42p cancels the correction coming from the 
partition function ([4T|) . Therefore, 



E (^H^m>r. 



(43) 



w=±i 



The first term in the r.h.s of Eq. has been eval- 
uated in Sec. IIII Al As it always happens in instanton 
physics^^, the derivative ipw{t) coincides with a zero 
mode of the fluctuation Sipw (t) . It is worthwhile to men- 
tion that only zero modes of fluctuations around instan- 
ton solution contribute to the non-pcrturbative renormal- 
ization of the polarization operator. The corresponding 
contribution is as follows (see Appendix C for details): 



W=±l 



In- 



T UT / (44) 

X e-^/^cos2nq + 0{ujl). 

From Eqs and (gi]) we obtain 

U{iujn) 2g2 Ec 

— ^' in — cos27r(7 



(J2 



^) 



(45) 



where constant D = (tt^/S) exp(— 7 — 1). 

The average charge on the island can be expressed via 
the partition function as 



T ainZ 
2Ec dq 



(46) 



Using Eq. (|¥T|) we find the following temperature and 
gate voltage dependence of the average charge in the one- 
instanton approximation: 



Q = q 



In — sin 2iTq. 



(47) 



Performing standard analytic continuation in Eq. (|45|) . 
we obtain the retarded polarization operator H^(a;) in 
the form of Eq. UT]) with 7ro(r) satisfying Eq. (gOl) and 



^,{T)^2^C'(-y-^~Dge 



■9{T)/2 



COS 2TTq 



(48) 



Finally, the average energy dissipation rate will be given 
by Eq. with function 



AT) 



9{T) 



(1 - i:).g2(r)e-ff(^)/2 cos 2719) 



(49) 
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In deriving this result we changed g to g(T) in the factor 
in front of the exponent in the r.h.s. of Eq. It is al- 

lowed by the accuracy wc arc working within. Result ((49|) 
asks for an interpretation. As expected, Coulomb block- 
ade manifests itself as a periodic dependance of dissi- 
pation A{T) on gate charge q. If we ascribe this de- 
pendance to the quantum resistance only, i.e., we write 
AiT) = C^RqiT) with Rq{T) following from Eq. (glD 
we face a paradox. It is believed that Coulomb blockade 
should suppress the tunneling of electrons between the 
island and the lead stronger for integer values of q than 
for the half-integer ones. Therefore, it would be natural 
to expect that Rq{T) is smaller at a half-integer value 
of q than at an integer one. The discussion above sug- 
gests that we have to conceive some kind of temperature 
renormalization of the gate capacitance Cg. 



D. Physical observables and gate capacitance 
renormalization 



As shown in Ref . [2^ , the proper physical observables 
for the Coulomb blockade problem are 



g'{T) = 47rlm 



dK^{ui) 



q'(T) = + Re 



duo 



j=0 



duj 



(50) 



j=0 



where the average charge Q is given by Eq. (|46p . the 
retarded correlation function K-'^{uj) is obtained from the 
Matsubara correlator 

K{t,2) = -|«(ri2)(e»[^(^^)"^(^^)l) (51) 

by standard analytic continuation. The physical observ- 
ables g'{T) and q'{T) describe a response of the system to 
a change in the boundary condition (|26p . One-instanton 
contribution to the physical observables g' and q' reveals 
their periodic dependence on the external charge q as^i^ 

g'{T) = g{T) (l - Dg{T)e-<^^'^^^^ cos 2TTq) , 

(52) 

q'{T) = q - —g^T)e-s<'^^^^ Sin27rq. 

Some remarks on the physical meaning of these quan- 
tities are in order here. In the perturbative regime g'(T) 
coincides with the rcnormalized coupling constant g{T) 
while q'{T) doesn't undergo any renormalization and 
coincides with the external charge: q'{T) = q. Thus, 
roughly speaking, we can think of them as the physi- 
cal observables corresponding to the action parameters 
g and q. The physics behind quantities (j50p becomes 
even more pronounced if we turn from a SEB to a sin- 
gle electron transistor (see Fig. [2]). In the absence of 
DC-voltage between left and right leads, a SET is de- 
scribed by the very same AES action (221)- (EZl) in which 
the bare coupling constant g = gi+gr- Here, gi/r denotes 



the dimcnsionlcss conductances of the left/right tunnel- 
ing junction. The quantity g'(T) then coincides with the 
SET conductance's*^ up to a temperature independent 
factor: 



G(T) 



91 gr 



h {gi + gr) 



r9'{T). 



(53) 



Expression for q' (T) is possible to write in terms of anti- 
symmetrized electron current-current correlatori ^^i'^^ 



q'iT) =Q-i 



25i5r dVdc 



(54) 



where Vdc denotes the DC-voltage between the left and 
the right leads and /(<) = dfid{t)/dt - the current opera- 
tor for the SET. 

For reasons to be explained shortly, it is natural to 
define the rcnormalized gate capacitance 



dq'iT) 
dUo ■ 



(55) 



According to Eq. ([5D)) . the quantity Cg{T) is different 
from the effective capacitance dQ/dUo which has been 
considered in the literature so far. On the perturba- 
tive level Cg{T) coincides with Cg; only instanton effects 
make it temperature and gate voltage dependent: 

Cg{T) ^Cg(l- 1.92 (r)e-3(^)/2 COS 2^<z). (56) 

Now, as usual, we plug in the bare capacitance Cg ex- 
pressed via Cg{T) into Eq. Wc see that instanton 
corrections cancel each other and the result ((49l) for the 
function A which determines the energy dissipation rate 
becomes 



A{T) 



27rC^^(T) 
9{T) • 



(57) 



With the same level accuracy wc can substitute 9'{T) for 
g(T) and obtain finally the following expressions for the 
energy dissipation rate and the admittance in the quasi- 
static regime 



= luj'C'g{T)Rq{T)\U^\^ Rq{T) 



eV(T)' 



(58) 
(59) 



giu,) = -^oJ^^+^Cl{T)Rq{T)u^\ 



Several remarks arc in order here. The results (j58p and 
([55)1 are valid in the weak-coupling regime: g'(T) ^ 1, in 
which the quantities dQ/dUo, g'{T) and Cg{T) are given 
by Eqs (gT]), ^ and respectively. Relations ^ 
and ([59)) fully describe the quasi-static dynamics of SEB. 
The energy dissipation rate factorizes into the product 
of well-defined physical observables in complete analogy 
with classical expression ([T]). The admittance behavior 
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..Ech{n,q) 




FIG. 3: (Color online) Charging energy Ech = Ec{n — qY as 
a function of gate charge q. 



is different from what we were expecting to get. Indeed, 
its imaginary and real components involve two different 
capacitances: effective capacitance dQ/dU^ and renor- 
malized gate capacitance Cg{T). Moreover, the temper- 
ature independent factor C/Cg survives in the real part 
oig[w). 



IV. STRONG COUPLING REGIME, g < f 

As follows from Eqs (|50p . the physical observables 
g'{T) and q'{T) are defined for arbitrary values of g. 
Therefore, it is of great interest to compute energy dis- 
sipation rate and the SEB admittance in the opposite 
regime, of small dimensionless tunneling conductance 
g ^ 1. The question we ask is whether the results 
(PI and dUD with the proper Cg{T) and Rq{T) hold? 
We mention that the case g <C 1 is a strong coupling 
regime from the field-theoretical point of view. In what 
follows, we compute energy dissipation rate by means of 
two different approaches. The first one is a refined field- 
theoretical method centered around Matveev's projective 
Hamiltoniaufi^ The second one is more straightforward 
approach of rate equations on which the 'orthodox the- 
ory' of Coulomb blockade was basedj^ We demonstrate 
how these two approaches beautifully complement each 
other. 



A. Preliminaries 

We center our effort around the most interesting case: 
the vicinity of a degeneracy point: q = k + 1/2 where k 
is an integer. Following Ref. [l^, the hamiltonian 
([6]) can be simplified by truncating the Hilbert space of 
electrons on the island to two charging states : with Q = 
k and Q = k + 1. The projected hamiltonian then takes 
a form of 2 x 2 matrix acting in the space of these two 
charging states. Denoting the deviation of the external 
charge from the degeneracy point by A: q = k + 1/2 — 
A/{2Ec) we write the projected hamiltonian as)^ 



H = Ho + Ht + AS^ 



where Hq is given by Eq. (|4|) and 

k,a 



(61) 



Here, S'^, = S'^ ± iS^ are ordinary (iso)spin 1/2 op- 
erators. The presence of small AC component in the 
gate voltage changes the parameter A according to: A 
A — {eCg/C)Uuj cos Ltit. This time the response of the 
system to AC gate voltage is determined by the isospin 
correlation function nf-(tj) (see Appendix A) which Mat- 
subara counterpart is given by 



n,(T) = (%s^t)s^o)). 



(62) 



The energy dissipation rate and SEB admittance can be 
expressed as follows 



(J2 

= ^c.Imnf(c.)|C/„P, 



-^c.^^f•(c.). 



(63) 



Therefore, we need to proceed with the computation of 

n.(T). 

To deal with spin operators it is convenient to use 
Abrikosov's pseudo-fermion techniquci^ We introduce 
two-component pseudo-fermion operators '0^, ipa such 
that 



(64) 



Pseudo-fermions bring in the redundant unphysical 
states when J2a V'iV'a > 1. To exclude these states one 
adds an additional chemical potential 77 to the hamilto- 
nian. It is necessary to set jy —00 at the end of any 
calculation. The physical partition function Z and corre- 
lators (C) can be found from the pseudo-fermionic ones 
as 



hm . 



(O) = hm 



d 



{0)pf 



(65) 



The elegance of pseudo-fermion technique lies in the fact 
that diagrams with pseudo-fermion loops vanish when 
one sets: r/ —00. 

Next, we plug representation ((64)) into the hamiltonian 
([SO)) , switch to Matsubara basis and integrate out elec- 
trons in the lead and the island. Done in the parametric 
regime ([TT]) this leads to the following effective action: 



drip 



cr.A 



dri dT2 a ( Ti 2 ) [ V' ( Ti ) cr _ -0 ( r 1 ) ] [ V' ( T2 ) cr+ V' ( T2 ) ] 



/3A^ 

iEc 4 



(66) 



FIG. 



Feynman rules for pseudo-fermion action; 



Here, ct^ stand for Pauli matrices and a± = {ax ±iay)/2. 
Action similar to Eq. ([66]) has been first analyzed by 
Larkin and Melnikov in Ref. [40| . In modern terminology, 
Eq. corresponds to the XY case of the Bose-Kondo 
model for the spin l/24ii^i^ Effective action is 
very suitable for our purpose since it's coupling constant 
(7 <C 1 justifying perturbative expansion. 

First, we establish the relation between pseudo-fermion 
and physical partition function. From Eq. (|65p . we find 



Z = lim ZpfC 



(67) 



Here, we denoted Gct(t) = —{Tr'ijja{T)-iija{0)) the exact 
pseudo-fermion Green's function. The Feynman rules for 
action (j66p are shown in Fig. 3) In the zeroth order in 
we obtain 



1 



l£r. 



Z ~ 2 cosh 



/3A 



2p/ = 1, (68) 



where e„ = 7rT(2n -I- 1) and ~ — r/ + crA/2. Spin- 
spin correlation function (|62p written in terms of pseudo- 
fermions becomes 
1 



n.,p/(r) = -(r.[^(r)a^^(T)][^(0)a^V(0)]), 



(69) 



where the average is taken with respect to action 
The physical correlation function is obtained from 
IVs^pf{iLUn) according to Eq. ([55]) . 



ni(^) e + uj 




FIG. 5: Feynman diagrams defining the polarization operator 
in the lowest order. 



Here, ^^{x) denotes the Euler digamma-function. The 
analytical continuation should be made with some care. 
We want to recover the retarded polarization operator 
n^(a;) which is regular in the upper half-plane of w. Since 
il'{x) has poles at a;„ = —n with natural n, the operator 
([70)) has poles in both halves of a complex plane. We get 
rid of superfluous ones with the help of identity: 

■0(z) Z)^ -TT cot TTZ. (72) 

Using a well-known relation for real x 



1 

2x 



■ coth 'KX 



(73) 



together with Eqs ([55)) and ([S5)) . we arrive at the follow- 
ing expression for the imaginary part of the polarization 
operator: 



imnf(.).|-jA ^.eoth: 

^ (T = ±l 



2T 



(74) 



A ■r-^ A — auj 
I coth — - — > coth 
2T ^ 

cr=±l 



2T 



, A 
tanh — -. 
2T 



B. Spin-spin correlaion function Ilf^^f{uj). First 
order in g. 



We start by calculating the polarization operator ([69)1 
in the lowest possible order of perturbation theory. 
It happened that the first non-trivial contribution to 
Iis^pf(iujn) came from the first order perturbation theory. 
The relevant Feynman diagrams are depicted in Fig. |5| 
The computation of Hs.pf{i^n) is rather straightforward 
and is presented in the appendix. The result is 



47r^ 



4±4?Z!^e^''sinh^, 

[iLOnY 2 



(70) 

where F^iu) is a regular in the upper half-plane of lo 
function: 



E 

a=±l 



(A- 



2'KTi 



icrAx 
2^) 



(71) 



Expression has a striking feature. It's divergent in 
the limit w — > 0. Indeed, 



imnf(^) 



9 



(3A 



inu) sinh /3A 



0. 



(75) 



The explanation is as follows. In essence the correla- 
tor ([62)) describes the noise of a fluctuating charge inside 
a metallic island. It was computed firstly in Ref. [il] . 
The author of Ref. [3| however obtained a different (reg- 
ular at a; — > 0) expression. He used a special type of 
analytical continuation which yields a symmetric noise 
{{nd{t) , fid{0)}) ■ We, on the other hand, are interested 
in its antisymmetric counterpart which is the response 
function (|16p . It is exactly this retarded antisymmetric 
function which is obtained via standard analytical con- 
tinuation procedure. 

The unphysical divergency ([75p comes from the non- 
trivial and essentially non-perturbative infrared structure 
of a polarization operator In what follows we 



prove that the partial summation of some infinite classes 
of diagrams resolves this singularity yielding the result: 



w + e 



Imnf(c^) 



(76) 



where z ~ gA at T = 0. As seen from Eq. ([75]) the limits 
Lj Q and g — > do not commute which explains how 
the artificial divergency in Eq. (|75p arises. Now we pro- 
ceed with a more accurate computation of the correlator 

nf(^). 



C. One-loop structure of the pseudo-fermion 
theory 

Throughout all our computation we will need some 
knowledge of the one-loop logarithmic structure of the 
pseudo-fcrmion theory. The bare Green's function is 
modified by the self-energy: 



Ga{ie„) = - 



1 



(77) 



The leading logarithmic approximation corresponds to 
one-loop renormalization. As it is known^, one extracts 
self-energy I]o-(j£n) from a self-consistent Dyson equa- 
tion: 

_3_ 

An' 



-j-T^\u.m\G-„{ien + iLOr,^). (78) 



Here, we introduce lo„^ = 2-KTm. The vital observation 
40ll43l43t is that the action can be renormalized with 



only one scaling parameter Z. Performing standard an- 
alytic continuation, we find^ldS 



Gf\e) 



Z{\) 



-1/2 . . E,, 



(79) 



max{T,|A|,|e|}- 



Here, = -?74-crA/2. g = gZ'^{\) and A AZ2(A) are 
the renormalized coupling constant and gap, respectively. 
As we see the energy Ec plays the role of a reference 
energy scale. The important feature of Green's function 
(|79l) is its acquired width 



cosh ^ 



r.(e) = ^(£-e-.)- . , 

sinh cosh ^ 



(80) 



According to the hierarchy of energy scales considered in 
the paper, Ec^ T and, therefore, the logarithmic correc- 
tions ~ glnEc/T are not small and require special care. 
To get read of large logarithms we change the reference 
scale of the field theory ([55)1 from Ec to T. With the 
help of result (|79p we may rewrite the theory in terms 
of renormalized fields and running coupling constants: 
ipa \J Z{X)tp^, 9^9 and A ^ A. The action then 
becomes 



(81) 



n(a,) = ^ 




exact vertex 



£, cr exact Green's 



function 

FIG. 6; Dyson equation for polarization operator Y\a^pf{iujn). 



where SSc.t. stands for the action of counter terms. It is 
responsible for a consistent regularization of higher or- 
der (in g) corrections to the physical obscrvablcs. Ac- 
tion ([8T|) is very suitable for our purpose. All large 
logarithms are absorbed into coupling constants and 
fermionic fields. This allows us to drop counter-terms 
in what follows. To bind the observables defined at the 
reference scale E,. with the renormalized ones we shall 
need to establish scaling of the pseudo-fcrmion density 
Ppf = '^^{ipa'^'a-) £md ^-component of the total spin den- 
sity s^pf = {l/2)J2a'^{M<y)- As follows from Rcf. 
the pseudo-fermion density p has no scaling dimension of 
its own: 



Ppf 



(82) 



where now the average is taken with respect to ac- 
tion ([5T|) . The total spin density Sz has the same struc- 
ture as the term proportional to A in action (j66p . There- 
fore, it must have the same scaling dimension: 



(83) 



where again the average is taken with respect to the ac- 
tion (j81|) . For completeness we present the rigorous proof 
of Eq. (j83p via Callan-Symanzik equation in Appendix D. 



D. Dyson equation for the spin-spin correlation 
function nf*pj(aj). 

The graphical representation of Dyson equation for the 
spin-spin correlation function 



Hs,p/(iw„) = —'^TaiiSk,iek+iLOn,ii^n) (84) 
X Gaiiek)Ga{iek + iujn) 



is illustrated in Fig. [HI Here, T^{iek,iSk + iw„,iw„) de- 
notes the exact vertex function. Performing the analytic 
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continuation in the spirit of Ref. [4^, we find (see Ap- pendix D for details) 



^ e + uj , £ 

tanli — — tanli — - 

2T 2T 



(85) 



+r«««(£, s + Lo)G^ie)G^{s + uj) tanh ^ - T^^^{s, s + uj, uj)G^{e)G^{e + lo) tanh i±^}. 



The most important task is to single out singular at 
Lu ^ and terms in ((85|) . We shall treat them 

separately to avoid divergencies. Firstly, we recall that 
we are interested in the quasi-static limit. Therefore, we 
shall proceed under assumptions to <^ max{A, T} but 
u ~ ^max{A,T}. It is intuitively clear that a singular 
contribution always comes from the first term in the r.h.s. 
of Eq. ([55]) which involves the product G^G^ . Indeed, 
we observe that the pole structure of G^G^ always leads 
to a singular denominator of the type (w + 2igT^) as a 
result of integration. This happens due to the proximity 
of poles in G^ and G^. In contrast, the other terms 
with G^G^ and G^G^ are regular at ,g = a; = and, 
therefore, free of divergencies. We may compute their 
contribution setting g ~ and safely expanding the re- 
sult in Lu. The integrand in Eq. ([SS)) also has a series of 
Matsubara-type poles due to the presence of hyperbolic 
functions. These poles lead to logarithmically divergent 



sums. The latter are controlled by the renormalization 
scheme. In our case all leading logarithms are absent. 
They have already been absorbed into renormalized con- 
stants g, and A by the proper choice of reference energy 
scale. Thus we can drop all divergent sums over Matsub- 
ara frequencies. 



Perfoming integration over e in Eq. (|85p and expanding 
in UJ where it is allowed, we are able to write down a much 
simpler expression for Il^pj{uj): 



E 



/3 



16 cosh^ I; 



2T 



1- 



AFtR 



- UJ,LlJ) 



u; + 2igr,i^^) 

(86) 

Now we need to compute the vertex function 
r^-^-^(e,e -I- uj,uj). The vertex function T„{iek,i£k + 
iuj„,iujn) satisfies Dyson equation 



J 



Ta{iek,iSk +iuJn,iuJn) = 1 + \^m\G-a{iek + iuJm)G-^a{i£k +il^m + iUJn)T-^a{i£k +iuJm,i£k +il^m +iUJn,iLOn) 



(87) 



r 



which is shown in Fig. [71 The details of analytical con- 
tinuation are described in Appendix D where we prove 
that the Dyson equation for the vertex r^^^ {e,e + Lu,Lu) 
becomes 



^ARR 



(e,e + uj,Lu) = 1 



9 



dx , 



^ARR. 



Stt / 27r 



(x, X + LO^Uj) 



xG'l,{x)Gl,{x + u:){x ~ e) 

X + UJ 



2 coth ■ 



2T 



tanh ■ 



2T 



tanh ■ 



2T 



(88) 



To solve it we have to make some self-consistent guess. 
The apparent difficulty is that apart from singular factor 
G^G^ the integrand in the r.h.s. of Eq. ((55)) may have 
unknown poles coming from the vertex function T^^^'. 
We conjecture that these poles result in the contribution 
of order of unity and are small comparing to singular 
contribution from the product G^G^. Hence we may 



perform an integral in the r.h.s. of Eq. 
at the following result 



^ARR 



{e,e + uj,uj) ^l- 



9 ^-a -£ 



and arrive 



cosh ^ 



sinh ^ 



2T 



cosh^ 
cosh ■ 



- { cosh 

- UJ, uj) 



2T 



2T 



-zc^ + 2gr„,(^ 



(89) 



We see that the solution does have an additional series of 
poles in the e plane. However these are Matsubara-type 
poles and are irrelevant as was argued above. Setting 
the external energy e = ^g. we obtain the self-consistent 
equation on r^^^(fcr,CCT +uj,uj). The solution reads 



J.ARR,, , . 1 [uJ + 2^-g{T.,-T,)\[uJ + 2^-gT, 



(90) 
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FIG. 7: Dyson equation for the vertex To-. 

Here, Fa = rrr{£,a) is the width of the Green's function 
defined in Eq. ((50)1 . 

CoUecting Eqs ((86|) and ([89|l we write down the resuh 
for the spin-spin correlation function 



9 



A 



47rTsinh^ 



, A 

— coth — 
27r 2r 



(91) 



FinaUy, taking into account Eqs ([67]) and ([68]) , we obtain 
the following result for the spin-spin correlator: 



nf(^) 



47r Tsinh^ 



5A , A 



(92) 



Here we restored factor Z** which provides a correct scal- 
ing dimension of spin fields according to Eq. ((83)) . 



E. The admittance and the energy dissipation rate 



With the help of Eq. ((55)) we obtain the admittance of 
the SEB for frequencies u; <C max{A,T}: 



A 



C 47r Tsinh 



Y + ^ coth ^ 



(93) 



The average charge Q and the physical observables g'{T) 
and q'{T) can be found from the pseudo-fermion the- 
ory ()66p if one substitutes the transverse spin-spin corre- 
lation function 

A'.(ri2) = -|a(ri2)(^+(n)5-(r2)) (94) 

for K{ti2) in Eq. ((SD))'^'^ In the leading logarithmic 
approximation, the average charge and the physical ob- 
servable g' are given bj^ 



Qm = fc+^-^tanhA, 



9'{T) 



-9 A 
2TsinhA 



(95) 
(96) 



The temperature dependence of the other physical ob- 
servable q' is as follows^l 



q'{T) = k+--- tanh — . 
^ ' 2 2 2T 



(97) 



To get the energy dissipation rate we expand expres- 
sion ((M)) in uj. Using the identity dA = —Z'^dU^/C and 
Eqs ((95)) -((97 )) . we obtain the energy dissipation rate and 
the admittance of the SET in the quasi-static regime: 

W^ = \u^Cl{T)R,{T)\U^\\ i?,(r) = -5-^, (98) 

g{u) = + ^Cl{T)R,{T)u\ (99) 

Here, the renormalized gate capacitance and the effective 
capacitance becomes 



CAT) 



dq' 



E, 



' 2 Tcosh^A' 



dUo ' 2 Tcosh^ A 



(100) 
(101) 



2T 

Several remarks are in order here. The results ((98)) 
and ((99)) are valid in the strong-coupling regime, g ^ 1 
and near the degeneracy point |A| ^ Ec- The accuracy 
we are working with (the leading logarithmic approxima- 
tion) allows us to make the following key observation. 
The expressions for energy dissipation rate ((55)) and ad- 
mittance ((M)) cast in terms of the quantities dQ/dUo, 
g'{T) and Cg{T) coincide with the ones obtained in the 
weak-coupling regime. It makes us suggest that the re- 
sults ((95| and ((99)) are valid for all temperature range 
Ec T ^ 6 and all values of g. 

We mention that formula ((55)) for Rq{T) is a truly 
non-perturbative in g result. Despite obvious complica- 
tions we overcame to obtain it, the expression for g'{T) 
(stripped of all logarithmic scaling) is the same as ob- 
tained in a much simpler approach of sequential tunnel- 
ing. This approach known as the 'orthodox theory' of a 
Coulomb blockade will help us to shed light on the physi- 
cal meaning of results (|98p and ((M)) . Further calculations 
are formulated in the language of rate equations^ which 
lie in the basis of the orthodox theory. 



F. Rate equations approach 

The rate approach is less general since it is essentially 
a Fermi golden rule approximation. It overlooks virtual 
processes and is unable to reproduce logarithmic scaling 
of physical observables. On the other hand rate equations 
are much easier to solve than corresponding Dyson equa- 
tions used above in a field-theoretical treatment. We are 
going to demonstrate that rate equations allow us to find 
the admittance for frequencies which are not restricted 
by the condition lo <C max{A,T} imposed by the field 
approach. Eventually we will conceive a prescription on 
how the admittance formula ()93p can be generalized for 
arbitrary (but still not very large w <C i?c) frequencies. 

As above, only two charging states of the island are 
counted. We denote them as follows: state corresponds 
to the average charge Q — k and state 1 corresponds to 
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Q = k+1. Probabilities for each state are denoted as po 
and pi which satisfy Po + Pi = 1- Master equation has 
the standard forni)^ 



Po 



-TioPo + Toipi. 



(102) 



Here, Fqi and Fio are tunneUng rates from and to the 
metalhc island, respectively. We should keep in mind 
that tunneling rates Tq^/iq are proportional to dimen- 
sional conductance of the tunneling junction g which is 
the expansion parameter of our problem. The average 
current through the contact is / = —po- Since we are 
interested in the linear response of the current to the 
time-dependent gate voltage U{t) = Uq + Uuicosujt, we 
expand the tunneling rates to the first order in amplitude 



01/10 



2C 



7oi/io('^)e" 



(103) 



-7oi/io(-^^)e" 



Then it is easy to find the following relation for the ad- 
mittance: 



c, 7io(a^)rE]i-7oi(u.)r?o 



The equilibrium tunneling rates are well-known^ 
r°win - 4^fcothA±i 



- 01/10 



Ait V 



2T 



(104) 



(105) 



We mention that up to the logarithmic corrections 
-'^01/10 ~ 2(7r±. A straightforward calculation of the tun- 
neling rates yields (see Appendix E) 



701/10 (t^) 



.9_ 



1 ± —F'^iuj) 

TTLO 



(106) 



where function i<fl(a;) was introduced in equation (|7ip . 
Plugging Eqs (fT051) - (fT06ll into Eq. (fT04|) wc arrive at the 
general expression for admittance: 



-icjcoth ^ - ^F^{uj) 



C47rcoth^ -iuj- 



0_ COtll nrri 



(107) 



In order to relate this result to field-thcorctical re- 
sult ((M]) . we expand the function F^{lo) to the first order 
in u): 



which is valid for u> <C max{A,T} and nearly repeats 
result for the admittance. The only difference is in 
the scaling factor Z which is absent in the rate equations 
approach. Now we may easily guess a prescription on 
how to generalize Eq. (|93|) for arbitrary oj. A correctly 
defined observable admittance ought to scale as Z^. It 
should also be expressed in terms of the renormalized 
parameters g and A only. This leads us to the following 
result 



C 47rcoth^ 



-iuj - 



-^'l^^ (110) 
coth 



which as we believe describes the admittance for lo ^ Ec 
in the strong-coupling regime g <C 1. Here, the function 
F^{lli) is given by Fji{uj) in which A is substituted for 
A. Finally, we mention that at finite frequency uj the 
parameter A which determines the scaling parameter Z 
in Eq. ([79)) should be modified as follows: 



A = In 



Er. 



max{r, |A|, 



(111) 



V. DISCUSSIONS AND CONCLUSIONS 

As we demonstrated in the previous sections the energy 
dissipation rate W^j is given by Eq. ^ with Rq{T) = 
h/e^g'{T) and Cg(r) dq'{T)/dUo in both weak and 
strong coupling regimes. We emphasize that the physical 
observables g' and q' are defined in terms of the correla- 
tion function of the phase field (^(r) of the AES model 
(see Eq. ([50| '). Therefore, they can be found, in gen- 
eral, not only in weak and strong coupling regimes but 
for arbitrary values of g and q. Hence, it is natural to 
assume that Eq. ^ as well as Eq. hold in general 
for a SEB under conditions of applicability of the AES 
model which are i?Th ^ Ec ^ T ^ (5max{l, g} and 
g/Nch ^ 1- Although, at present we are not able to 
prove this conjecture we believe strongly it is true. 

Originally,— the physical quantity q'{T) has been in- 
troduced for a SET and its physical meaning was inter- 
preted in terms of the average charge on the island and 
the antisymmetrized current-current correlation function 
(see Eq. (j54p . If we introduce non-symmetrized current 
noise in the SETiLi^ 



Si{uj,Vdc) 



dte.- 



(112) 



F^{lo) = ITTLU 



2rsinh 



, A ' 

„ ■ — coth — 



(108) 



Plugging this into p04p we get the familiar expression: 
Cg g PA —iuj 



g{u) 



C 47r sinh /3A -i^ + |^ coth A 



(109) 



then wc can present Eq. (|54p as 



q = 



{9l+9rf 

2Tr gigr 



-p.v. 



duj dSiiio, Vdc) 



UJ dVd, 



Vdc=0 



(113) 



Therefore, to measure q'{T) two separate experiments 
are needed: the measurement of the average charge on 
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the island at Vdc = and the measurement of the non- 
symmetrized current noise 5*7(0-', V^c)- Although exper- 
imental designs probing the non-symmetrized current 
noise have already been proposed^ and measurements 
have recently been taken from a number of electronic 
quantum devices^, it is still a challenge. Our present 
results indicate that the quantity g' can be related to the 
rcnormalized gate capacitance Cg{T). Namely, Cg{T) — 
dq' /dUo, provided that result ^ holds in general (not 
in a weak and strong coupling regimes only). The latter 
capacitance can be extracted from measurements of the 
energy dissipation rate and the SET conductance. 

Recently, the energy dissipation rate of SEB has 
been addressed experimentally via radio- frequency reflec- 
tometry measurements (by sending a continuous radio- 
frequency signal to the device).— The temperature and 
external charge dependences of the quantity [uj^A{T)]~-^ 
were studied. The latter quantity was refered to as the 
'Sisyphus' resistance by the authors of Ref. [l3|- In the 
experiment the tunneling conductance was estimated to 
be equal g w 0.5 such that the SEB was in the strong 
coupling regime. In Ref. the 'Sisyphus' resistance 
was estimated theoreticall y w ithin the rate equation ap- 
proach (see Eq.(4) in Ref. [i2|)- Their result corresponds 
to our result (|109p for the admittance. However, our fi- 
nal result for the admittance (|110[) is more general then 
Eq. (4) of Ref. [l3| ■ The latter does not take into account 
not only the logarithmic renormalizations of the SEB pa- 
rameters but also deviation of the function F^'{u!) from 
the linear one. Although the values of the SEB parame- 
ters reported in Ref. are such that the difference of 
the scaling factor Z from unity is several per cent, log- 
arithmic renormalizations in the expression for the ad- 
mittance yield noticable effect. This is shown in Fig. [S] 
In addition, the function F^{iLj) can be written as the 
linear one only for frequencies lo <^ max{A, T} which is 
not the case for the low temperature data of Ref. p^ . 
Therefore, the experimental data of Ref. [l3| needs to be 
reanalyzed with the help of Eq. (|110p . 

The authors of Ref. [l^l claim that their results for the 
'Sysiphus' resistance indicate the violation of the Kirch- 
hoffs laws. They argue that the admittance they measure 
does not correspond to the equivalent circuit of SEB with 
bare values of the gate capacitance Cg and the tunneling 
conductance g. However, by the same logic one could 
claim the violation of the Kirchhoff's laws in measure- 
ments of the SET conductance G{T) because it is dif- 
ferent from gigr/{gi + gr)- Our results imply that the 
energy dissipation rate (inverse of the 'Sysiphus' resis- 
tance) in the SET can be obtained from the Kirchhoff's 
laws if one substitute Cg and g for Cg{T) and g'{T) in 
the equivalent circuit. As one can see from Fig. [HI the 
energy dissipation rate is maximal for A = which cor- 
responds to the half-integer values of the external charge 
q. It occurs because the larger value of A is, the larger 
the ratio rpj^/F^Q becomes. We remind that F[jj^y-^(, is the 
transition rate from (to) the state with Q = fc + 1 to 
(from) the state with Q = k. The increase of A makes 
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FIG. 8: The dissipative part of admittance of the SEB at 
fixed cj as a function of A. Three plots using three different 
formulae are presented. Gr is given by pU7p . Go is given 
by (fT09| . Gf is given by (nO\i . We use g = 0.5, Ec = lOT, 
and Lu = O.ST. 

the system less probable to be excited in the state with 
Q = k + lhy the time-dependent gate voltage and, there- 
fore, reduces energy dissipation. Of course, this physical 
explanation is strongly based on the model of only two 
charging states involved. It is valid at 5 ^ 1. How- 
ever, at g'{T) ^ 1 the energy dissipation rate has the 
maximum at half-integer values of the external charge q 
as well (see Eq. (03) )• This result cannot be explained 
by arguments based on the 'orthodox' theory since there 
arc no well-defined charging states in the weak-coupling 
regime. 

The dissipation caused by the electron tunneling is not 
the only one that occurs in the set-up. Intrinsic elec- 
tron transitions inside the metallic island cause an addi- 
tional internal energy loss. This mechanism is, in fact, 
the origin of metallic conductivity. This sort of dissipa- 
tion ought to be mainly a classical effect. It corresponds 
to the radiation of energy by a metallic particle placed 
in the quasi-stationary electric field. The classical dissi- 
pation can be conveniently characterized by two limiting 
regimes: the low-frequency ohmic loss and high frequency 
non-ohmic radiation (skin-effect): 

5te^ gta^n 

^3/2 (114) 

gte/ \u}oJ 

Here, a — e'^/hc is a fine structure constant, gte^/h is 
an internal (Thoulcss) conductance of the island, R - its 
characteristic size and luq is the separating frequency. To 
elucidate the parametric conditions under which c^uan- 
tum dissipation Wui due to presence of the tunneling 
junction dominates over the classical one we make nec- 
essary estimates. Quantum dissipation can also be split 
into ohmic and non-ohmic limiting regimes. The cor- 
responding separating frequency is denoted as fi. We 
are concerned with simple estimates only and drop weak 
log-corrections in all formulae for the quantum case. The 
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FIG. 9: Schematic diagram for comparison of quantum and 
classical mechanisms of the energy dissipation. The quantum 
dissipation dominates in the filled region. 



TABLE I: Description of functions for Fig. |9] 











LOO 
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^2 






ii 

$2 


1 (^2^ 

,,7/2 J 


7,7/2 yS2i J 



results are most transparently explained via phase dia- 
gram which is presented in FiglS] supplemented by Ta- 
bles U and |TT1 

In the fully coherent case, the admittance of SEB was 
studied in Ref. 0| by means of the S'-matrix formalism. 
It was shown that the SEB admittance can be presented 
in accordance with its classical appearance but the 
definition of physical quantities comprising it becomes 
different. In Ref. it was derived that the gate ca- 
pacitance Cg and the tunneling resistance R should be 
substituted by the mesoscopic capacitance and the 
charge relaxation resistance respectively. However, 
according to our results, although being applicable in the 
fully incoherent case, the SEB admittance in the quasi- 
stationary regime involves two capacitances: the effective 
capacitance dQ/dUo which controls the imaginary part 
of Q{uj) and the renormalized capacitance Cg{T) which 
together with Rq(T) determines the temperature behav- 
ior of ReC/(w). It is the effective capacitance that corre- 



TABLE II: Description of parameters for Fig. |9] 





n g 


> 1 


gE,/h 




g < 1, A < T 


gT/h 


.g < 1, A > T 


gA/h 


AT A/T 



sponds to the mesoscopic capacitance C^. The appear- 
ance of the effective capacitance dQ/dUo in the imagi- 
nary part of the admittance is dictated by conservation 
of charge via the Ward identity . We expect that the 
SEB admittance should involve two physically different 
capacitances in general. Recently, the SEB admittance 
was studied with the help of the S'-matrix formalism in 
the incoherent case also In particular, it was predicted 
that in the fully incoherent regime and at low tempera- 
tures the charge-relaxation resistance Rq = h/{ge'^). It 
is at odds with our result that Rq = h/{e'^g'[T)) since 
at low temperatures g'{T) can be very different from g 
(see Eqs. ([52]) and ([96]) ). The reason behind this dis- 
crepancy is as follows. Coulomb interaction in Ref. [l3| 
was accounted for on the level of classical equations of 
motion only, which was the conservation of charge. In 
the mean time quantum fluctuations of charge are signif- 
icant throughout all our analysis and there is no obvious 
justification to take them negligible. 

To summarize, we have studied the energy dissipa- 
tion in a single electron box due to a slowly oscillating 
gate voltage. We focused on the regime of not very low 
temperatures when electron coherence can be neglected 
but quantum fluctuations of charge are strong due to 
Coulomb interaction. We considered cases of weak and 
strong coupling. In both cases we found that the energy 
dissipation rate is given by the same expression involv- 
ing two physical observables g'{T) and Cg{T). Our result 
for the energy dissipation rate can be obtained from the 
SEB equivalent circuit if one substitutes g' {T) and Cg (T) 
for g and Cg, respectively. We strongly believe that the 
universal expression we found for the energy dissipation 
rate is valid for an arbitrary value of the tunneling con- 
ductance. 
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VI. APPENDIX 
A. Energy derivative 

Here we relate dissipation in the system to various field 
correlators in weak and strong coupling regime. 

a. Weak Coupling, g ^ \ We want to express cor- 
relator p6)) through AES effective phase (/)(t). We want 
to be rigorous and introduce Kcldysh contour. We split 
all the fields into upper and lower components (±) in 
correspondence to the halves of the Kcldysh contour. 
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FIG. 10: Keldysh contour. 



The action of the system is split as well S = S-^ — S- 
and the partition function of the system reads: Z = 
J ^(y9±e*^['^±l = 1. The average electron density is found 
as: 

(A.l) 

C / 6S ■ ^ ' 



Here, we introduced classical and quantum components 
for bosonic fields: 

Ua,c^\{Ug+ + Ug^), U,^, = \{U,+ -Ug^), (A.2) 

and Ug{t) ~ Uq + U^cosut. To get rid of quartic 
Coulomb terms we introduce Hubbard-Stratonovich bo- 
son fields: V- on each part of the contour and make 
fcrmion gauge transformation: 

rfa. -da.e-'/o^"'^*. (A.3) 

The transformed terms ^o, Sc, St take the form: 

/oo 

+ E / diJidt-e^^Ad^^dt, 

^ i nOQ nOO 
Sca^^ dt + Cg V^UgAt) dt, 

^ J-oo J -co 

a.k 

We see that the source term Ug,a enter Sea only, hence 
we regroup it in a more suitable form: 

Sc ~ Sc+ Sc— = 

/oo />oo 
VcVg dt + V2Cg / {VcUg,q + VgU g ^ dt . 

-OO J —oo 

(A.5) 

Here, Vc,q = "^(^+ ± VL). Next we find the physical 
electron density from (jA.l^ : 



(A.4) 



C 



(A.6) 



V2' 



We then expand e''^ to linear order in physical field Uc 
The result reads 



c. 



2 poo 



/oo 
nn{t-t')Ug{t')dt', (A.7) 
-OO 



where Hfl(i - t') = iC^ {Vcit)Vq{t')) . Coupled with ^ 
it gives (pg)) . 

Using (jA.6[) we also write down the formula for the 
effective capacitance dQ/dUo of the SEB: 



dQ__ . H^,(0) 
dUo~ ' C ■ 



(A., 



b. Strong coupling, g <C 1 We proceed in complete 
analogy with the previous case. Using Hamiltonian (|60p 
we obtain 



dH 



Co 



dUgl C 



(S^t)). 



(A.9) 



Keldysh technique gives 



(S'it)) 



c 



i{S',{t)Sl{t'))U,{t')dt' . (A.IO) 



Introducing spin-correlation function 

Hf (t) = i{smsim 



(A.ll) 



we recover dissipation expression (|63[) with spin correla- 
tor n^(ci;) playing the role of polarization operator. 



Admittance 



The admittance is defined as: 

Q{uj)e 



SUg{t') 



2tt 



(B.l) 



We introduce the tunneling current operator using hamil- 
tonian (1210: 

/ = i[H, ^ 4d„] = i J2 tka^aldo. + h.c. (B.2) 



To find the average current we insert the necessary source 
term into the action: 



1 /"OO 
Ss = - 

2 



I{t)K{t)dt, = 



i dK{t) 



(B.3) 



While taking a functional integral along Keldysh contour 
we keep a quantum component of K{t) field only. We 
make the usual rotation in the fermion basis: 



1 



1p± = -^{■4j2±'4'i)- 



(B.4) 
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Here ip = {a^, da)'^ . After the rotation and gauge trans- 
formation (jA.3[) the source and tunneling terms take the 
form: 



St + Ss = I dtilj^[T^s{^) + -J^5{f)ji's 



Jq Jc 
Jc Jq J 



(B.5) 



Here, indices c and q denote classical and quantum com- 
ponent of a corresponding physical value, i.e. Jc,q = 
\{J+ ± J~) and Ao-, are matrices in a island-lead 
space: 



tkae~ 

Ak^'^' 

itkaS 



If a 



(B.6) 







It's possible to get rid of highly non- linear source term 
(|B.3p by a suitable change of field variables. Indeed, one 
can easily check that up to linear order in k: 

(B.7) 

The same property holds for all the elements of matrices 
Tjs, J-yS- By making a change: 



we put the whole k -dependence into gaussian part of the 
action. Then 

Ss = - j «(0(^V^c + CgUc)dt. (B.9) 



The average current (|B.3p reads 



(B.IO) 



Using (jA.Sp and to linear order in Uc{t) wc find the cur- 
rent to be 

(/) = CgU, + ^{(P,) + iCCg J U,{t')dt'{(p,{t)^,{t')). 

(B.ll) 



The admittance becomes 
Hence, 

C RQg{uj) 



ImHj^ (w) 



Cg UJ 



(B.12) 



(B.13) 



In the case of spin variables (strong coupling) we can eas- 
ily get the analogue of formula (|B.12[) for the admittance 
using the same steps. This way we establish the relation 
between admittance and spin polarization operator H^ 
quoted in the main body: 



c, 



c 



(B.14) 



where n^(w) is given by (jA.lip 



C. Instanton contributions 

1. Massive fluctuations 

We expand the fluctuating field 5if{T) in the basis of 
eigenfunctions: 5ip{t) ~ '}2m^m'~Pm{'T), where the basis 
reads {u ^ e^'^''^^)^ 



1 ~ UZ 

\ 1 — uz 

ip-rn[T,z) = -—^ , 



m > 2 
m > 2: 



(C.l) 



u — z 
u 



Here, <y5±i(T, z) are field zero modes. Then the correlator 
reads 



(B.8) {%5^{t)5^{t')) =tY, J S!z^-m{T,z)^^{T',z) 



2n 



1 



1 



gt^m-i g{m~l)T' 
d^z ■ ■ Z_ 



TO > 



\z\ < 1 



Here, Vi/Vq is the ratio of fluctuation determinants. 
Some care should be taken when regularizing them. We 
used the scheme proposed iit^^ 



Vq 27r3T' 
After simple algebra wc obtain 



(C.2) 



21n(l - s) 




+ s 
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Expanding this expression into Teylor series over s we 
get 



I = ^ns"ln^, 



n=l 
oo 



n— 1 /c— 1 n— 1 

oo 



III = 2^s", 



n=l 

oo n 



ri=l fc=l 



Here, i7„ is harmonic number. The contribution of gaus- 
sian fluctuations into the correlator becomes 



_L-e"'-"'"'{%Sv(rmr'))„ 



n=l 



n=l 



S 



(C.3) 



Now we make analytical continuation of Fourier- 
components into the region n <^ I. We are interested 
in linear in n term. 



Extracting linear part and summing instanton and anti- 
instanton terms we obtain 



■ in — I cos znq 



4ttT T 



(C.4) 



which does cancel partition function renormalization 

(ED. 



2. Zero modes 

The corresponding single instanton configuration 
reads: 



Lpw = 2'kTW{- 



— z 1 — zu 
The correlator is given by 



W = ±l. (C.5) 



^{|z2|V + |zYV"}. 



(C.6) 



The corresponding Fourier-component is as follows 

V r U^l'"' 

{TMrMr'))„ ^ e~'^/^87r^Tcos27rq-^ / ^-^dh. 

T^o J I-fI 

(C.7) 

Expanding it in n ^ 1 to linear order we reproduce (|44p . 

D. Computation of polarization operator 

1. The lowest order 

First we notice that l{—uj„) = II(cl;„). Thus we will 
drop any odd function of tOn while calculating l{LUn). The 
analytical expression for diagram I (see Fig[5]) reads 



1 



1 



i{ek + LOn + fi,n) - i-a iSk - S.a ' 

Performing the sum over fermion frequencies we get 

l(Wnj 



47r 



(iw„)2 



1 



Act + i{LL!n + flm) Act + f2„ 
1 



Act + i{u)n + ^m) (Act -|- Z^m)^ 



where nf{x) = l/(e''^ + l) is Fermi distribution function. 
Simple algebra shows that I(a;„) -|- II(ti;„) = — III(lj„). 
Thus, taking the limit 77 — > —00 we obtain 



I(w„) + II(w„) + III(c^„) 



^e*'sinhA^|n„ 



(A + if7„, + iuJn){A + iVL^Y 



Now it's clear that the sum over Vim can be taken in terms 
of digamma functions. The answer is given by Eq.(|70p. 



2. Callan-Symanzik equation for (s^) 

The anomalous dimension 7 of operator s^j is intro- 
duced as 



Z^s;;(A,5) = 4^(A,5,A) 



(D.l) 



where Z is given by ([75]) and A is a cut-off: A ^ Ec- To 
extract 7 we write down the corresponding CS-equation 
for the Green's function: Fpf{A,g,A) — 5 X^cr '''(V^o■V'cr)■ 
The tree-level Fpf reads 



F,y(A)=-e^"sinh— . 



(D.2) 
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Following general RG-philosophy and with the help of 
(|D.1[) we write the corresponding CS-equation for func- 
tion F(A, g, A) in the form: 



(D.3) 



where the corresponding /3 - functions are given by 



Pa 



2n' 



27r2 



(D.4) 



The term with Pg always contains extra g and can be 
dropped in the leading order. Using action (|66|) we work 
out the last term: 



dlnZ g 



(D.5) 



d In A Att'^ ' 
To find 7 we need to get F in the next to (|D.2[) order: 



A 



F.fiA, g, A) = -e^" sinh i± (l - ^ In ^ 
+ e''.^eoshAin^. 



(D.6) 



Here, e - is a characteristic scale of interaction. Plugging 
(ID3)) and (|DJl) into (|D3|) we find: 



7 = 2. 



(D.7) 



a retarded vertex function LOn > 0. Next we define three 
vertex functions in accordance with the structure of cuts: 

r^'^^^(z, z + iu, iuj) if Imz > 0, 

r^^^(z,z + iw,icj) if -iw„<lmz<0, (D.IO) 

V^^^^^z^z + iuj^iuj) if Imz < -icjn. 

The general expression for lla{iLOn) then becomes: 



-E 

4 ^ 



ds £ 

— — tanh— ro.(£,£ + iuJn,i^n)Ga{e + zcj„)Go-(£). 

Q lOTTl ZJ 

(D.ll) 

The contour C is shown in FigHT] As usual the integral 

© 

yrrr I 

Im(e) -_0_ ' 




lm{£ + iuj) — 0\_ 




3. Exact expression for polarization operator 

In order to work out the polarization-operator diagram 
in FiglS] we follow the scheme proposed by Eliashberg 
[45| . First we establish the analytical properties of ver- 
tex function T{z, z + iuJm i^n) as a function of complex 
variable z. The operator expression for the vertex func- 
tion reads 

r,(Ti - T,T2 - r) = {%i}„{t)^a{t)i^a{tl)^a{t2))- (D.8) 

Its Lehman representation is as follows 



T„{z,z+ioj,iuj) = T^ J2 wr„ 



'ikln 



Z + LOkl 



Z + iuj — UJir, 



Wjm -iuj \ Z + UJ„rn. 

Wil^k = {l\i^a\n){m\i;,\k). 



Z + iuj - UJln 



(D.9) 



Complex calculus teaches us that the sum (|D.9p defines 
a function with two horizontal cuts: Im(z -I- iuj) = and 
Im(z) = 0. For simplicity let's restrict our attention to 



FIG. 11: Contour for polarization operator YI{lu). 

over large circle vanishes and we are left with integrals 
over different branches: 



1 



de , e 

tanh — 

47ri 2T 



^AAR. 



(e - iLUn,s,iujn)G^{e)G^{e - iw„)|. 



Making analytical continuation icun ^ Lo+iO we get result 
183. 



4- Dyson equation for the vertex 

Following the same scheme as in the previous section 
we derive the expression for the vertex function. The 
contour G depends on the type of the vertex we need 
to get from (jST]) . The contour for the vertex r^^-'^ is 
depicted in Fig[T2l The result reads 



^ARR 



{e,e + uj,uj) = l-I-fll-III, 



(D.13) 
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Im(z + ie) = 




Im 2 = I L 

^ 1 «' 



As usual regularization scheme allows us to drop these 
sums. The integrand of term I however contains G^G^. 
As a consequence it is singular at g — s- as explained 
in the main body. This way we recover Eq. 



Ivciiz + iu + ie) = 0\ 




FIG. 12: Contour for the vertex function r-*-^-^. 



where 



1 = j ^ ^GUx)G\{x + c.)r:4^^(^, x + u^u:) 



2ilmafl(x — e) coth 



X ^ £ 

2T 



. X X + LU 

- tanh —aR{x - e) + tanh - e) 



II 



dx 



ranix ~ e)G^^{x)G'l^{x + uj) 



X r 



III 



RRR 



(x, X + uj,uj) tanh 



2T' 



dx 



UAix ~ e)G^Jx)G^^ix + Lu) 



X r 



AAR 



{x, X + Ll},Ll}) tanh 



2T 



(D.14) 



Here, function a(z) is an interaction propagator whose 
Matsubara counterpart is shown in FigH) As usual it 
has a cut Imz = which allows to define two functions: 



_IUJ _IUJ 



(D.15) 



The integrands entering terms II and III are explicitly 
analytical in the upper and lower halves of the complex 
plane respectively. Consequently we may turn the corre- 
sponding integrals in to sums over Matsubara frequencies 
i£„. Next one can easily prove the following identities: 



^RRR 



^AAR 



^ARR 



■^ARR 



{ie„ - uj,ien,uj). 



(D.16) 



This way we drastically simplify our Dyson equation by 
rewriting it entirely in terms of a single vertex T^^^. 
Next, 



{iSn - S,-a){i£n + - £,-cr) ' 

47r 



(D.17) 



lEn — ijj — e 



E. Rate probabilities 

To work out rates r° and 7 we follow standard scheme. 
We introduce Heisenberg ^/j-operators according to 



(E.l) 



Then the matrix elements in the basis of filling numbers 
become 

{0\H^)^Y.^O\d^\^)e-''^'-+^'>K (E.2) 



In an ordinary fashion we change the Hamiltonian by 
gauge transformation of fermion fields (path-integral ap- 
proach is implied): 



(E.3) 



where U{t) = Uuj cos ujt. Now the whole U (i)-dependance 
is transferred into the tunneling part of the Hamiltonian: 

Ht^Y. tkcald^e'^ f ^(*)''* + h.c. (E.4) 



Let us compute rate rio(t). The initial and final states 
read 



K) = |fc,iv), 

\f)=diak\k,N). 



(E.5) 



Here, k (N) is the number of electrons in the island (lead). 
As usual, S'-matrix formalism gives the necessary ampli- 
tude in the form: 



(E.6) 



Now we substitute / U{i)dt = (U^/uj) sinwt and tunnel- 
ing Hamiltonian assumes the form: 



Aio{t) = -t{f \ / Ht{t)dt\i) 

J — 00 

= -i(z|4d„ / Ht{t)dt\i). 



Ht = y^^tkaa\da{l 



^^^sinc^t) +h.c. (E.7) 



The detailed-balance relations for probability rates read 

(E.8) 



rOi(A)==r?o(-A), 

7oi(<,A) = -7io(t,-A). 
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Plugging (jE.7p into (|E.6p and integrating one gets the 
amplitude of transition — )■ 1 



Aio{t) ^ tl,{l - n^)nke'^'--"'+^^' 
\ 1 



£fc - - A + iO 2Cw 



Ek-ek' - A-u; + iO Sk - Sk' - A + uj + iO 

(E.9) 

Squaring it, taking thermal average and integrating we 
get the full expression for the probability in the linear 
response regime: 



Wio{t) 



9 



„2At 



2/5" - 1 (S - A)2 + A2 



Cuj 



1 



(E.IO) 



where g is defined in (|12p . Now wc find the transition 
rate as a derivative of a transition probability Tio{t) — 
dWio{t)/dt and the following expression for 7s: 



710 (w) 



"27r 



ds 



2m e^^" - 1 ■ 



1 



1 



s-A + iOs-A-w-iO 

1 1 
s-A-iOs-A + cj + iO 



(E.ll) 



The integrand converges very well in the complex plane 
and the integral can be easily taken: 



710 (t^) = - 

00 



2ti 



1 



A + uj 



1 



ILOr, 



. ( ILUri — I 



A — LO iljJn — A + UJ 



(E.12) 



Expressing the sum in terms of digamma functions we 
get (fTOe]) . With the help of ((RT21) and (jKS]) one can 
establish the following useful identity: 



710 (tj) - 701 (t^) = — . 



(E.13) 
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